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Abstract 


Simple  shearing  deformations  are  analyzed  as  a  means  of  in¬ 
vestigating  the  effects  of  strain  hardening,  thermal  softening,  and 
strain-rate  sensitivity  on  the  formation  of  shear  bands.  Linear 
stability  analysis  is  conducted  for  perturbations  from  homogeneous, 
time  varying  deformations.  Finite  difference  solutions  are  obtained 
for  the  fully  nonlinear  system  of  governing  equations.  Velocity 
boundary  conditions  are  emphasized  although  solutions  for  the  case 
of  stress  conditions  are  also  presented.  Three  types  of  constitu¬ 
tive  models  are  considered:  the  power  law  model,  the  Arrhenius 
law,  and  the  Bodner-Merzer  model.  Inertia  and  heat  conduction 
are  included.  Qualitative  features  of  plastic  response  for  the 
class  of  problems  considered  are  illustrated  by  eighteen  numerical 
solutions.  General  conclusions  on  the  critical  conditions  for  the 
stability  of  the  solutions  are  difficult  to  draw  in  view  of  the  many 
features  of  the  problem  that  influence  the  formation  of  shear 
bands.  Low  strain  hardening,  high  thermal  softening,  and  weak 
strain-rate  sensitivity  all  contribute  to  shear  strain  localization. 
The  size  of  an  initial  perturbation  is  important  in  the  stability 
of  solutions  for  the  case  of  velocity  boundary  conditions,  but  it 
is  less  important  for  stress  boundary  conditions  because  of  the  in¬ 
herent  inhomogeneity  of  dynamic  solutions  in  the  latter  case. 
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I .  INTRODUCTION 


Localization  of  plastic  deformation  in  the  form  of  shear  bands 
is  observed  frequently  during  tensile,  compressive  and  shearing  defor¬ 
mation  tests.  Shear  bands  occur  in  a  wide  range  of  applications  such 
as  ballistic  impact,  dynamic  fracture,  cryogenic  behavior  of  materials, 
high  velocity  forming,  machining  and  grinding. 

Current  understanding  of  the  phenomenon  is  not  fully  established; 
however  a  common  view  is  that  at  high  strain  rates  the  occurence  of  shear 
bands  is  related  to  thermal  softening  features  of  material  behavior.  That 
is,  shear  bands  are  thought  to  initiate  under  conditions  of  nearly 
adiabatic  deformation  when  thermal  softening  overcomes  strain-harden¬ 
ing  so  that  the  rate  of  change  of  the  flow  stress  becomes  equal  to 
zero  [1].  At  this  stage  of  deformation  the  flow  stress  decreases 
as  temperature  increases,  resulting  in  inhomogeneous  deformation. 

This  inhomogeneous  deformation  results  because  the  heat  generated  by 
plastic  working  is  greatest  in  regions  of  highest  strain  rate.  This 
leads  to  further  softening  of  the  material  in  these  regions,  and  the 
cycle  is  repeated  until  the  deformation  is  localized  in  the  form  of 
shear  bands.  The  existence  of  material  or  geometric  inhomogenieties 
appears  to  be  necessary  for  triggering  this  mechanism.  The  exact 
nature  and  size  of  such  inhomogeneities  are  not  fully  understood; 
however,  it  is  believed  that  these  defects  play  a  major  role  in  shear 
band  formation. 

Early  observations  and  discussions  of  shear  bands  are  attributed 
to  Zener  and  Holloman  [2] ,  who  observed  a  white  etchinq  band 
in  an  0.25  C,  low  alloy  steel  subjected  to  a  rapid  punching  test. 

They  showed  that  a  local  temperature  of  1000°C  is  reached  with  a  total 
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shear  strain  of  (5.0) .  A  large  number  of  similar  experiments  are  cited 
in  a  survey  by  Rogers  [3].  Strain  rate  sensitive  materials  are  viewed 
as  being  less  sensitive  to  localization  [4] .  Some  metals  show  very 
stable  behavior  in  their  superplastic  ranges  where  they  are  well 
characterized  as  nearly  Newtonian  fluids.  Rate  independent  models  of 
plastic  response,  without  thermal  softening,  exhibit  shear  strain 
localization  when  the  effects  of  geometric  softening  due  to  lattice 
rotations  and/or  non-normality  of  the  plastic  strain-rate  to  the 
yield  surface  are  included.  These  effects  have  been  presented  fully 
in  a  recent  review  by  Asaro  [51 . 


Experimental  investigations  of  shear  bands  often  are  concerned 
with  metallurgical  effects  on  the  onset  of  localization  [6] ,  and  with 
measuring  the  strain,  strain-rate  and  temperature  at  which  bands  are 

formed  [7], [12].  Analytical  discussions  of  shear  bands  seek,  for 
assumed  constitutive  models,  the  critical  conditions  which  character¬ 
ize  the  instabilities  reported  experimentally.  Most  of  the  assumed 
constitutive  models  are  based  on  a  macroscopic  point  of  view  due  to 
the  associated  complexity  of  plastic  deformation  at  the  microscopic 
level.  Moreover,  the  pre-localization  constitutive  equations  are 
assumed  to  hold  when  localization  occurs.  Several  authors  attempted 
to  obtain  a  simple  stability  criterion  for  adiabatic  plastic  deforma¬ 
tion  by  means  of  a  linear  perturbation  analysis  [8],  [9],  [10] ,  [ll] . 
However,  the  use  of  constant  coefficients,  evaluated  from  the  homo¬ 
geneous  solution  at  the  time  the  perturbation  occurred,  imposes  a 
severe  restriction  on  the  validity  of  such  results.  In  other  analyses, 
computer  simulations  have  been  presented  for  the  fully  non-linear 
system  [12 J ,  rl 3] . 
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This  work  presents  a  general  analysis  of  shear-strain  local¬ 
ization  in  simple  shearing  deformations  of  thermal-visco/plastic 
materials.  The  system  of  governing  equations  for  a  one-dimensional 
simple-shearing  deformation  is  formulated  in  Chapter  II.  An  exact 
solution  for  quasi-static  deformation  with  no  heat  conduction  is  pre¬ 
sented  in  Chapter  III  for  materials  with  temperature-dependent 
viscosity  and  no  strain-hardening.  In  Chapter  IV  a  linear  stability 
analysis  is  presented  in  which  the  time-dependence  of  the  coefficients 
in  the  linearized  equations  is  retained.  A  stability  criterion  moti¬ 
vated  by  a  bifurcation  type  of  instability  is  also  presented. 

In  Chapter  V,  a  numerical  simulation  of  the  fully  nonlinear 
system  of  equations  is  presented  for  a  number  of  constitutive 
equations.  These  constitutive  equations  exhibit  strain-rate 
sensitivity,  strain-hardening  (softening)  and  thermal  softening 
(hardening)  that  are  representative  of  the  response  of  metals.  The 
major  objective  underlying  such  numerical  simulations  is  to  study  the 
material  response  to  small  initial  perturbations  of  the  homogeneous 
solution.  This  work  serves  as  a  means  to  obtain  a  better  insight  into 
the  phenomenon  of  shear-strain  localization,  and  to  reconsider  some  of 
the  common  ideas  concerning  shear  band  formation. 


II. 


PROBLEM  FORMULATION 


Consider  a  simple  shearing  deformation  of  an  infinite  plate  at 

constant  strain  rate.  The  plate  has  thickness  H  in  the  y-direction 

as  shown  in  Fig.  II. 1.  The  lower  edge  (y  =  O)  is  fixed,  while  the 

upper  edge  (y  =  H)  has  a  velocity  (v  =  v  )  in  the  x-direction. 

o 

Assume  that  all  physical  quantities  are  uniform  along  the  x  and  z 
directions  so  that  the  deformation  depends  only  on  one  space  coordin¬ 
ate.  Under  these  conditions,  the  only  non-zero  displacement  component 
is  in  the  x-direction,  and  the  only  non- vanishing  stress  components 

are  a  ,  o  ,  a  ,  a  .  The  equations  of  momentum  balance  reduce  to: 
xy  xx  yy  zz 
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where  0 ^  and  c^  are  the  components  of  the  Cauchy  stress  tensor,  v 
is  the  particle  velocity  in  the  x-direction;  and  pQ  is  the  mass 
density.  Hereafter  the  shear  stress  will  be  denoted  by  t.  The 
total  strain  is  assumed  to  be  composed  of  an  elastic  plus  a  visco-plastic 
part,  i.e. 

Y  =  Y®  +  YP  (II. 2) 


where  the  elastic  part  of  this  decomposition  is  given  by 
e  T 


Y 
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(II. 3) 


in  which  G  is  the  elastic  shear  modulus.  Differentiating  (II. 2)  and 
(II. 3)  w.r.t.  time  and  replacing  the  total  strain  rate  by  the  velocity 


gradient  one  obtains 
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(II. 4) 


The  heat  Q  generated  by  plastic  deformation  is  assumed  to  be 
related  to  the  plastic  work  by 


52  =  6t  i£ 

3t  p  3t  Pl  3t 


(II. 5) 


where  the  factor  3  is  taken  to  be  approximately  0.9.  Heat  transfer 
in  the  material  is  assumed  to  satisfy  Fourier's  law  of  heat  conduction 


h  =  -  k 
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30 

3y 


(II. 6) 


where  h  is  the  heat  flux  in  the  y-direction,  k  is  the  thermal 
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conductivity  and  9  is  the  temperature.  From  (II. 5)  and  (II. 6)  the 
energy  balance  equation  can  be  written  as 


PoC  ft  =  k  —  +  6t 

°  3t  3y2 


*r 


3t 


(II. 7) 


where  c  is  the  specific  heat  per  unit  mass. 

The  remaining  equation  is  the  constitutive  equation  that 
characterizes  the  material  response  to  the  given  loading.  We  consider 
materials  which  exhibit  strain-rate  sensitivity,  strain  hardening, 
and  thermal  softening  or  hardening.  Thus,  in  simple  shear  the  flow 
stress  is  assumed  to  be  given  by  a  relation  of  the  form 

,P 


t  =  tp  (0,y  ,  r ) 


( ;p= !x  ) 

at  1 


(II. 8) 


The  governing  system  of  equations  used  for  modelling  the 
features  of  shear  strain  localization  in  thermal  visco-plastic  materials 
can  be  summarized  as  follows: 
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T  = 

ip(0»  rp. 
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(11.12) 
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(II. 12a) 

The  complete  initial-boundary  value  problem  is  defined  once  initial 
and  boundary  conditions  are  prescribed.  These  conditions  will  be 
given  when  numerical  solutions  are  discussed  (chapter  v) • 

Dimensionless  Form 

It  is  useful  to  rewrite  the  system  of  governing  equations  (II. 9, 
12)  in  a  dimensionless  form.  The  dimensional  quantities  in  the 
system  (11.9,12)  will  be  relabeled  by  a  superposed  ~  ,  and  then  the 
unlabeled  quantities  will  denote  the  corresponding  dimensionless 
quantities. 

The  dimensionless  quantities  are  defined  by 

t  =  t  ,  y  =  y/H  ,  x  =  t/t  , 


k  0  t 
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°  T  H2 
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where  the  subscript  "o"  refers  to  quantities  evaluated  for  the 
homogeneous  solution  at  the  time  a  perturbation  is  introduced,  while 
the  superscript  "o"  refers  to  evaluation  for  the  homogeneous  solution. 

A  *“  D 

The  characteristic  time  t  =  1/y  is  the  time  required  to  obtain  a 

o  o 

unit  shear  strain  at  the  strain  rate  y  P  .  The  coefficients  V,S  and 

o 

T  are  defined  by 


V  =  1±  s  =  T  =  li 

V  3t  '  S  3yP  30 

They  are  measures,  respectively,  of  the  strain  rate  sensitivity, 
strain  hardening  (S  <  0) ,  and  thermal  softening  (or  hardening) . 

The  dimensionless  quantities  defined  earlier  are  chosen  such  that 
once  introduced  in  the  system  of  governing  equations  (II. 9. 12)  ,  the 
resulting  system  will  be  exactly  the  same  except  that  all  quantities 
are  dimensionless.  The  dimensionless  form  will  be  used  in  the  remain¬ 


der  of  this  work. 


III.  ANALYTICAL  SOLUTION  FOR  A  SPECIAL  CASE 


Analytical  solutions  for  the  nonlinear  system  (II. 9-1 I. 12)  are 
difficult  to  obtain  due  to  nonlinearities  associated  with  the 
constitutive  equation  (11.12).  Examples  of  the  nonlinear  constitutive 
equation  are  the  enqpirical  relation. 


Qv  ,  pm,  3y  ,n 
T  -  Uo  6  (l  )  (  ) 


(III.l) 


where  V*o»v,  m,  n  are  material  parameters,  and  the  Arrhenius  law 
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(III. 2) 


where  u)Q  is  the  frequency  of  vibration  of  dislocations  as  they  attempt 
to  overcome  obstacles,  t (yP)  is  the  stress  barrier  amplitude  -  which 

is  a  function  of  the  plastic  strain  in  case  of  a  strain  hardening 

* 

material;  w  is  the  Activation  volume;  is  Boltzmann's  Constant. 

For  physically  reasonable  models  these  quantities  satisfy. 


a)  >>  1  ,  i  <  x_(  y  )  , 

• 

b 


(III. 3) 


The  function  x  (y*)  is  assumed  to  be  positive,  monotonic  and  to  have 
13 

a  saturation  limit  x  °°  as  y^  -*■  °°  ;  furthermore,  its  derivatives  are 
assumed  to  vanish  as  y^  ->  »  . 

In  order  to  be  able  to  obtain  an  exact  solution  that  exhibits 
some  features  of  the  solutions  of  (II.9-II. 12),  consider  the  following 


special  case: 


quasi-static  deformation,  i.e. 


at 


o 


•  no  heat  conduction,  i.e.  k  =  O 

•  no  elastic  effects,  i.e.  G  ->■  00 

•  constant  shear  stress  boundary  conditions,  i.e. 

t(  o  » t)  =  t(H,  t)  = 

•  Non-Newtonian  fluid  with  temperature  dependent  viscosity,  i.e. 


T 


v  3v  n 
P  (6)  (|^ 

o  ay 


n  >  o 


(III. 4) 


The  effects  of  thermal  -  softening  (hardening)  are  obtained  for 
negative  (positive)  values  of  the  power  \>.  Although  this  model 
does  not  represent  solids  of  engineering  interest,  it  is  attractive 
because  it  allows  an  analytical  solution  that  gives  insight  into 
instabilities  of  solutions  of  the  system  (II.9-II.12)  as  well  as 
providing  a  reference  for  checking  the  validity  of  the  solutions 
obtained  by  a  linear  perturbation  technique  and  by  a  finite  difference 
method.  Under  these  assumptions  the  system  of  equations  (II.9-II.12) 
reduces  to 


3x 

3y 

0 

(III. 5) 

9YP  _ 
3t 

3v 

3y 

(III. 6) 

36 

3t 

3v 

ri  T  3y 

(III. 7) 

T  = 

U  (9)V(|^)n 
o  3y 

(III. 8) 

where  r,  =  B/p  c*  From  (III. 5)  and  the  stress  boundary  conditions 
1  o 


the  shear  stress  is  given  by 
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T(y,t)  =  i. 


(III. 9) 


which  is  kept  constant  throughout  the  deformation.  Eliminating 
3  v  „ 

(— )  between  equations  (III. 7)  and  (III. 8)  one  obtains 
3y 


=  ne 

at  1 

,T1.  n  . 


(III. 10) 


n 

where  r\  =  r^x^( — )  is  a  positive  constant  and  v  =  —  .  Integration  of 

o 


equation  (III. 10)  gives 


0 (y ,t)  =  {f (y ,t,v)  } 


(III. 11) 


where  f(y,t,v)  =  0Q(y)V+^  +  n(vhl)t.  (III. 12) 

This  expression  for  the  tenperature  6(y,t)  can  be  used  to  obtain 
the  following  expressions  for  other  quantities  of  physical  interest 


If 


(y,t)  = 


yp  (y,t) 


n  |  f (y, t, v)  | 
p  ~  de  (y)  -I  r 

['.wv  '  { 

(jr)"'  {  }*' 


(III. 13) 

1<Mv) 

f (y , t, v) [  (III. 14) 


(III. 15) 


where 


(III. 16) 


Analysis  of  the  exact  solution: 

The  functions  f(y,t,v)  and  <$>(v)  characterize  all  the  features 
of  the  solution  to  the  problem.  The  function  <j>(v)  is  negative  for 
v  >  O  and  v  <  -1,  while  it  has  positive  values  for  -1  <  v  <  O. 

The  function  f(y,t,v)  is  always  positive  for  v  >  -1.  For  v<-l  there 


:v:  vS-s- 
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is  a  critical  time 


(6  (y )) 
_ o _ 

nj(v+i)| 


(III. 17) 


at  which  f  becomes  zero;  f  is  negative  for  t  >  t  Three  types  of 
behavior  can  be  distinguished. 

(i)  Thermal  hardening  for  v  >  0 

For  v  >  0,  <j>(v)  is  negative  and  f(y,t,v)  is  positive.  Then,  from 

(III.  13-III.  15)  the  quantities  f-7  ,  77-  ,  and  vanish  as 

at  a  t  ay 

t  -*■  «  .  Such  solutions  are  strongly  stable. 


(ii)  Thermal  softening  for  -1  <  v  <  0 


For  v  in  this  range  both  <j>(v)  and  f(y,t,v)  are  positive. 


The  quantities  ,  -  ,  —  ucuuuc  uuuuunucu  »  u  ^  w  . 

3t  3t  3y 

The  solution  grows  in  time  as  shown  in  Fig.  (Til. 2)  for  the  plastic 


30  3£_  39 


become  unbounded  as  t  +  ° 


strain  rate.  The  shape  of  the  plastic  strain  rate  distribution 
varies  slowly  with  time  at  small  times  and  becomes  independent  of 
time  as  t  -*■  00 .  There  is  no  tendency  for  localization  of  shear 
straining.  Based  on  the  uniformity  of  the  strain  distribution; 
solutions  such  as  those  shown  in  Fig.  dll. 2)  are  said  to  "grow" 
rather  than  to  become  unstable.  The  term  "instability"  is  reserved 
for  bifuration  phenomena  in  which  the  homogeneous  and  perturbed 
solutions  grow  at  different  rates  so  that  differences  between  these 
solutions  become  large  at  late  times  as  shown  in  Fig.  (III. 3). 

This  notion  of  instability  corresponds  to  experimental  observations 
of  shear  bands  with  localized  deformations  which  are  large  relative 
to  deformations  in  neighboring  regions. 


»■-*-  V  v-.'V  '■'l'  v-,  C.  v1  ^ 


(INITIAL  PERTURBED  SOLUTION) 


-  Uniform  Growth 


(iii)  Thermal  softening  for  v  <  -1 


In  this  regime  <j>(v)  is  negative;  f(y,t,v)  is  positive  for 

t  <  tcr  and  becomes  negative  for  t  >  t  .  All  quantities  in  (III. 11) 

and  (III . 13-III . 15)  become  unbounded  as  t  +  t  .  Solutions  exhibiting 

cr 

such  behavior  are  termed  unstable,  although  both  the  homogeneous  and 
perturbed  solutions  become  unbounded  at  the  same  time. 

Concluding  Remarks  on  Exact  Solutions 

1.  The  exact  solution  shows  the  possibility  of  a  type  of 
instability  in  which  the  solution  becomes  infinite  in  a  finite 
time. 

2.  Less  thermal  softening  is  required  to  obtain  unstable 
behavior  for  weakly  strain-rate  sensitive  materials  (i.e.  n  <<  1) 


since  Condition  (iii)  is  V  <  -n. 


IV.  LINEAR  STABILITY  ANALYSES 


In  view  of  the  extreme  difficulty  of  solving  the  non-linear 
problem  it  is  of  interest  to  conduct  stability  analyses  on  a  linear¬ 
ized  version  of  the  problem.  Several  authors  have  discussed  this 
approach  [ 9  ] ,  [11 ] ;  however,  for  most  of  this  discussion  the  homo¬ 
geneous  solution  which  is  to  be  perturbed  is  held  constant  at  the 
values  corresponding  to  the  instant  at  which  the  perturbation  is 
introduced.  This  assumption,  which  is  made  to  obtain  linear  equations 
with  constant  coefficients  severely  limits  the  applicability  of  the 
analyses.  Extensions  of  the  analysis  to  cases  of  time  -  varying 
coefficients  have  been  presented  by  Clifton  [14]  and  by  T.  J.  Burns 
[1C],  who  employed  the  direct  Lyapunov's  method  for  a  thermo- viscous 
material,  (i.e.  elasticity  and  strain  hardening  are  neglected) 

In  this  study  we  seek  a  rigorous  linear  stability  analysis  to  the 
general  problem  defined  by  Eqs.  (II.9-II. 12) .  Then,  we  present  a  modified 
stability  criterion  that  extends  the  validity  of  the  linear  approach. 


IV. 1.  Linearized  system  of  equations 

We  consider  the  linearization  of  equations  (II.  9-12)  about 
the  homogeneous  (parallel  flow)  solution.  We  seek  solutions  of  the 


form: 


t  (y,t) 

=  x°(t)  + 

6x (y ,t) 

v(y,t) 

ii 

< 

o 

>< 

r+ 

+  6v(y,t) 

YP  (y,t) 

=  yP°(t) 

+  6yP(y,t) 

(IV.  1) 

9  (y*t) 

=  0° (t) 

+  60 (y ,  t) 

where  6x ,  6v,  6yP,  60  are  the  differences  between  the  solution  x,  v, 
YP»  9,  that  includes  the  effects  of  a  perturbation,  and  the 
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O  O  p'J  o 

homogeneous  solution  x  ,  v  ,  y  ,  0  .  The  homogeneous  solution  is 
assumed  to  be  perturbed  at  (t  =  0)  by  a  small  periodic  fluctuation 
in  the  y-direction.  We  substitute  solutions  (IV. 1)  into  the  non¬ 
linear  system  (II.  9-12),  and  retain  only  first  order  terms  in 
6t,  6v,  6yP,  <50.  In  this  way  we  obtain  the  following  linear  system 
of  partial  differential  equations  in  6x,  6v,  6yP,  and  60. 


3  ,,  ,  1  3  . 

—  (6v)  =  —  —  ( 6  t  ) 

3t  p  3  y 

o 

h  -  F  l6vl  -  5  F 


(IV-2) 


(IV- 3) 


—  (60)  =  — -  (68)  +  r1x°  —  (6yP)  +  r1  (IV-4) 

3y 


—  -  \T°X.T  j. 


(6yF)  =  V  6x  +  T  60  +  S  6y 


° 


(IV- 5) 


where 


_  k 

o  ~  p  c 
o 

In  these  equations  the  coefficients  V°,  T°,  S°  are  defined  by 


v°  _  f  Jl_ 
I 

T°  =  Ji. 

30 


is  a  measure  of  strain-rate  sensitivity; 


T  =  [^"gg  j  is  a  measure  of  thermal  softening  (T  >0), 

or  hardening  (T  <  0)  -  thermal  softening  is 
the  usual  case  for  crystalline  solids; 

.  A  o 

O  3r 

S  =  — —  is  a  measure  of  strain-hardening,  (or  softening 

3yp 

if  negative). 

The  functions  V°,  T°  and  S°  are  evaluated  at  the  homogeneous  solution 
at  time  t.  hence, these  are  time  dependent  coefficients.  For  convenience 
of  notation,  we  suppress  the  "6"  in  the  following  analysis  (e-g-  6t  -*•  x)  . 
We  rewrite  equations  (IV. 2,5)  in  matrix  form  for  a  solution  vector 


18  - 


V  _■  _»  j»  .• 


The  resulting  system  of  equations  has  the  form 


(IV. 6) 


L  (w)  =  A  u>  +  A  u),  +  Ay  w,  +  A  y  w,  —  0  (IV. 7) 

-  t  ~  y  ~  yy 


where  a  comma  denotes  differentiation  with  respect  to  the  following 
subscripted  variable  and  the  coefficient  matrices  are  given  by: 


D  0  0  0 

o 


■r1  (T0V°+<f>°) 


o  o 

-V  T 


0.0 

-r1x  S 


t  0-00 
A  =  G 

0  0  10 


O  C  O  ] 


(IV. 8) 


O  O 


0  0 


0  0 


0  0 


O  0  -r 


The  coefficient  matrix  A,  is  time-dependent,  but  it  does  not  depend 
on  the  spatial  coordinate  y. 

We  introduce  the  Fourier  transform  defined  by 


F(€,t)  = 


e  F  (y , t)  dy 


(IV.  9.a  ) 


r\-‘ 

K‘* 

/  ,♦ 

l'  •  «  * 


sms* aaa 
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F(y,t)  = 


elCy  F(C,t)  dC 


Then  taking  the  Fourier  transform  of  L(u)) 

% 


8?  <«•« 


C(£,t)  iu(5,t) 

Of 


or,  suppressing  the  dependence  on  C, 

djij 

—  (t)  =  C(t)  u)  (t) 

dt  t 


O  we  obtain 

a, 


(IV. 9. b) 


(IV. 10) 


(IV. 11) 


where  C(t)  is  given  by: 


C(t)  = 


0 

-H/P 

o 

0 

0 

i£G 

-GV° 

-GT° 

-GS° 

0 

r^  (t°V°+4>0) 

°m°  r2 

V  T  “V 

rlT°s° 

0 

v° 

T° 

S° 

(IV. 12) 


Equation  (IV. 11)  represent  a  system  of  linear,  first  order  ordinary 

a  __ 

differential  equations  for  the  perturbation  vector  u)(£, t) . 

' \j 

IV. 2)  Solution  of  the  Linear  System: 


Analytic  solution  of  the  system  (IV. 11)  can  be  obtained  in  an 
iterative  way  as  follows: 


w (o)  +  [c(n)  to  ]  dn 

*  >o  ' 


=  u(0)  +  \C(T])  l  to (O)  +  [C(x)to]  dx] 


fn  - 

[C(x)to]  dx]  ^  d  n 


Continuing  the  iteration  process,  we  reach  the  solution: 

to  (t)  =  S(t)  w(0)  (IV. 13) 


where,  the  matrix  S (t)  is  known  as  the  matrizant  of  the  differential 
equation,  or  the  "solution  operator",  which  is  given  explicitly  by 


S(t)  =  I  + 


ft 

ft 

c(n)  dn  + 

c(n) 

o  J 

O  J 

C(x)dx  dn  + 


(IV. 14) 


Or,  in  compact  form, 

oo 

S(t)  =  £  C(t)(n) 

n=o 


(IV. 15) 


where 


C(t) 


(n) 


ft 


c(n)  c(n)(n_1)  dn 


and 


C(t) 


(o) 


I  . 


(IV. 16) 


(IV. 17) 


The  solution  expressed  by  (IV. 14, 17)  is  the  general  solution  of 
equation  (IV. 11);  however  in  the  special  case  where  C(t)  commutes 
with  C(t)^  the  solution  operator  S(t)  can  be  expressed  concisely 
in  the  form 


S(t)  =  exp  [ 


c (n)  dn  ] 


(IV. 18) 


(see  Appendix  I)  . 


Stability  of  the  Solution 


Consider  the  solution  in  the  general  form  (IV. 14)  which  repre¬ 
sents  an  infinite  series.  Taking  the  norm  of  each  side,  one  obtains: 


I  I  c(t) (n) ||  <  l  | |c(t) 

n=0  n=0 


However , 


'"'ll  -  Ilf- 


c(n)  c(n) (n-1)dn| 


and  for  finite  time  intervals  (0,t),  we  can  write 


C(n)  C(n) 


tn'1>  I  Idn  £  fl 


I c ( n)  I  |  |c (n)  (n-1)  dn 


|s(t)||  $  l  +  l  j  1 1 c ( n)  1 1  |  |c(n) 

n=l  I  J0 


(n-1) 


dn  \  .  (iv. 19) 


Let  us  assume  that  sup| |c(n) | |  exists  and  is  equal  to  a  finite 


value  p,  say.  In  this  case  we  obtain: 


|  S (t)  |  |  <  1  +  pt  +  ^ 2,  +  . . 


[S(t>  |  |  £  etsuPl|c<n>!l]t.  Vi)  e  (0>t] 


(IV. 20) 


From  (IV. 20)  we  conclude  that  the  series  S(t)  is  uniformly  and 
absolutely  convergent  for  all  values  of  t  as  long  as  the  maximum 
(supremum)  norm  of  the  coefficient  matrix  C (t)  is  bounded  for  any 
finite  time. 


IV. 4.  Stability  criterion  for  thermal-visco-plastic  deformation 

The  solution  (IV.13)  gives  a  direct  relation  between  the  vector 
of  initial  perturbations  and  their  corresponding  values  at  a  later 
time  (t  >  0).  The  argument  usually  used  to  determine  the  stability  of 
such  solutions  is  based  on  the  signs  of  eigenvalues  of  the  solution 
operator  S (t) .  Such  a  criterion  serves  to  determine  whether  the 
initial  disturbances  grow  (or  decay)  with  increasing  time.  However,  as 
mentioned  in  Chapter  III,  we  seek  a  bifurcation  process  for  which  the 
perturbed  solution  grows  much  faster  than  the  homogeneous  solution. 
Also,  we  require  the  unstable  growth  to  be  confined  to  some  region  of 
the  slab  thickness.  Motivated  by  this  notion  of  instability,  and  in 
order  to  distinguish  between  a  uniform  growth  Fig.  (III. 2)  and  an  un¬ 
stable  deformation  Fig.  (III. 3),  we  suggest  the  following  criterion. 

The  system  of  equations  (IV. 7)  is  called  "well-posed"  and  its 
solution  given  by  (IV.13)  is  called  stable  if  there  exist  constants 
k  and  a,  both  positive,  and  independent  of  initial  values  o>(0)  such 
that: 

sup  |  |  S  ( t )  |  |  <  ic  exp  (at).  (IV.  21) 

The  quantity  sup||s(t)||  in  (IV. 21) is  the  norm  | j S  (t) | |  defined  in 
(IV. 23) .* 


*If  (IV. 21)  holds  for  ||s(t) 

defined  by  j|s{t)||  =  max. 

x^O 


I  ,  then  it  also  holds  for  the  L~-norm 

l|s(t)x|l2 

-  for  each  t. 


m 


. • .  * . 
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The  notion  of  well-posedness  (IV. 21)  implies  continuous  dependence 
on  the  initial  data.  Moreover,  this  estimate  (IV. 21)  implies  that  a  stable 
solution  can  grow  exponentially  as  long  as  it  can  be  bounded  by  a  suitable 
choice  of  (<,a)  for  any  finite  time  (t)  and  finite  wave  number  (?) .  The 
inequality  (IV. 21)  suggests  much  weaker  conditions  for  stability  than  the 
requirement  that  all  the  eigenvalues  of  s(t)  satisfv  X.  <  0.  An  essential 

x  l  - 

feature  of  the  notion  of  instability  adopted  here  is  the  requirement  that 
there  are  solutions  that  become  unbounded  in  a  finite  time. 

In  order  to  apply  the  criterion  (IV. 21)  we  note  from  (IV.  20)  that 
(IV. 21)  is  satisfied  if  there  exists  an  a  >  0  such  that 


sup  I |c(n) | |  <  o.  (IV. 22) 

n€(0,t] 


That  is,  the  criterion  (IV. 21)  is  satisfied  if  the  norm  of  C(t)  is 


bounded  for  any  finite  time.  To  express  this  constraint  in  terms  of 

the  elements  of  C(t)  we  choose  sup  | |c(n) | |  to  be  defined 

n€[o,t] 

as  follows: 


n 

sup  j |c(n) | |  =  | | C (t) | | —  —  max  I  |c..|.  (IV. 23) 

n€[0,t]  l<i<n  j=l  13 


Applying  (IV. 23)  to  the  coefficient  matrix  C(t)  given  by  (IV.  12)  ,  we 
find 


=  sup  1 1 c  (n) 1 1 
n€[0,t] 


max(R1,R3,R4) 


where  the  row  sums  are 


R2  E  |G|[|V°|  +  I s° 1  +  I T° I ]  +  I GC | 

R3  =  I  rx  1  ‘  [  |t°V°  +  <j>°|  +  |t°S°|3  +  |r1T°T°-ro?2| 


i  | V° |  +  | S°  J  +  | T° | 


Since  for  constant  strain  rates  $  the  time  dependence  in  the  expressions 
for  the  row  sums  is  restricted  to  the  functions  T°,S°,V®  and  t°, 
we  can  write  the  stability  criterion  (IV. 22)  in  the  form 


{  | n)  1  +  |s°(n)|  +  |v°(n)|  }  i  a 


(IV. 24a) 


| x° (n) |  •  {  |T°(n)|  +  | s° (n) |  +  |v°(n) |  }  <_  6  (iv.24b) 

provided  that 

x° ( p)  >  ,  yn£[0#  t]  ,  (t  <  °°)  ,  a  >  o  ,  B  >  o  , 

O 

where  the  restriction  on  t  (n)  is  imposed  because  the  shear  stress  must 
be  greater  than  some  yield  value,  say  for  constitutive  equations  of 
the  form  (11.12)  to  apply. 

The  stability  criteria  (IV. 24)  have  been  derived  for  the  case  of  finite 
wave  numbers  £.  This  restriction  can  be  removed  by  considering  directlv 
the  rate  of  growth  of  the  I,  -norm  of  the  solution  ou  of  the  system 
(IV. 7)  .  (See  7g?pendix  IV) 


IV.  5.  APPLICATIONS  OF  THE  STABILITY  CRITERION 


IV. 5.1.  The  quasi-static  model  (III.5-III.8) 

The  analytical  solution  for  this  model  is  given  in  Chapter  (III) 
For  the  constitutive  law  (III. 8) 

o 

0 

S  *  0 

1-n 

„  0  x.  n 

V°  =  rf  =  (— )  (— )  (6°) 

Mon  yo 


-v 


Substituting  the  values  of  t0,T°,S°  and  V°  in  the  inequality  (IV. 24) 
we  obtain 


where 


(V  |v|Tl0°(t)-(V+1)  +  90(t)"V  <  a. 


(IV. 25) 


Substituting  the  homogeneous  solution 


0  (t)  =  [1  +  n(v+l)t] 


into  (IV. 25)  gives  the  stability  requirement 


(■”)  "  /  ^  +  |v|T1f(t;0)~1  <  a 


(IV. 26) 


where 


f (t ; v)  =  1  +  n(v+l)t. 


Obviously,  the  inequality  (IV. 26)  is  only  violated  for  v  <  -1.  In 
this  case  the  function  f  becomes  zero  at  a  critical  time 


fccr  n  |  \>+l  | 


(IV. 27) 


Therefore,  we  conclude  that  the  deformation  described  by  this  model  is 
stable  for  v  >  -1.  This  is  the  same  conclusion  obtained  from  the 
exact  solution  (III. ll-III . 16) .  Furthermore,  the  linear  stability 
analysis  predicts  the  same  type  of  singularity  as  that  obtained  for  the 
exact  solution,  as  well  as  the  same  critical  time.  Thus,  the  criterion 


(IV. 21)  appears  to  have  promise  for  characterizing  the  stability  of 


homogeneous  deformations. 


IV. 5. 2.  General  polynomial  constitutive  law  with 
velocity  boundary  conditions 

We  consider  the  general  system  of  equations  (II.9-II.12)  with  the 
constitutive  equation  (III.l)  written  in  the  form 

<f>  =  [7-  e"V(YP)"m]1/n  (IV.  28) 


A  homogeneous  solution  for  this  model  (neglecting  elasticity)  has 

been  given  by  Clifton,  et  al  [14] .  The  temperature  and  stress  in  this 

solution,  obtained  for  a  constant  strain  rate  <j>  ,  are 

0 


6°(t)  =  {yf  (  v,m)  •t^1+m) 


•r» 


1 

l-v 


(IV. 29) 


T°(t)  = 


0  v 
i®  it>] 


m 


(IV.  30) 


where , 


v  =  .  (1+m+n) 

y  -  v0*o  >  ° 


f(v,m)  = 


l-v 

1+m 


—  _  ,  n+m  _  _ 

w  =  vQtQ  >  0, 


For  the  constitutive  equation  (iv.28)  the  coefficients  T°,S°,V0 


are 


_,0  "0  .v. 

T  - (— ) 


(IV. 31a) 
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n  S0tJ 


(IV. 31b) 


12.  (i-) 

n  CT 

T 


(IV. 31c) 


From  (IV.29-IV.31)  the  stability  criterion  (IV. 24)  becomes 


~  •  U  +|m|2-  (6°)  Vtm"1+{v  |  y  (0°)  V’*1tm}  <  g  (IV.  32) 


From  (IV. 32)  the  stability  criterion  is  satisfied  as  long  as  neither 


0  v  0  v-1 

of  the  functions  (8  )  ,  (0  )  becomes  unbounded.  Such  singularities 


can  occur  only  for  f(v,m)  <0  and  at 


1/1+m 


y  j  f  (v,m) 


(IV. 33) 


Stability  boundaries  based  on  the  sign  of  f(v,m)  are  shown  in 
Fig.  (IV. 1). 

a 

The  function  f(v,m)  is  negative  in  regions  (I)  and  (III)  defined 
by 


I: 

V  >  1 

and 

m  >  -1 

III: 

v  <  1 

and 

m  <  -1 

Regions  I  and  IV  are  of  little  interest  in  applications  since  materials 
are  generally  thermally  soft.  Region  III  is  seldom,  if  ever, 
applicable  because  few  materials  exhibit  the  strong  strain  softening 
that  is  characteric  of  this  region.  Region  II  is  characteristic  of 
most  metals.  Thus,  it  appears  that  the  stability  criterion  is  satis¬ 
fied  for  most  metals  that  could  be  modeled  by  (iv.28).  In  view  of  the 
fact  that  Region  II  contains  points  corresponding  to  both  strain 
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softening  (m  <  0)  and  thermal  softening  (v  <  0) ,  it  appears  that 
stability  in  the  sense  of  (IV. 21)  is  satisfied  for  a  wide  class  of 
materials . 

Consider  the  term  in  {  }  in  equation  (IV. 32);  and  denote  it  p(t) , 
then  equation  (IV. 32)  becomes 

§ 

(— )  P(t)  <  a. 
n  — 

If  P(t)  is  bounded,  say  by  P  ,  then  the  solution  to  the  perturbed 

TuclX  • 

8t 

problem  grows  no  faster  than  eP  with  g  =  <f>n  P  „  .  Whether  such  a 

^  nictx  • 

n 

solution  to  the  perturbed  problem  grows  or  decays  depends  on  the  thermal 
softening  and  strain  hardening  parameters,  v  and  m,  respectively.  In 
either  case,  the  rate  of  growth  or  decay  is  affected  by  the  strain-rate 
sensitivity  parameter  (n) .  Small  values  of  n  (i.e.  weak  strain-rate 
sensitivity)  allow  high  rates  of  growth  or  decay. 

If  P(t)  is  unbounded,  then  the  problem  (IV. 11)  is  not  well-posed. 


IV. 5 . 3  Bodner-Merzer  Constitutive  Model 

Consider  the  system  of  equations  (II.9-II.12)  with  the  constitu¬ 
tive  model  (introduced  by  Bodner  and  used  by  Merzer  [13]) ; 

T  (Wp> 

T  =  t(yp,  e,  wp)  =  — -  [f (n) ]c  (IV. 34) 

/3 
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Neglect  elastic  strain  rates  in  (II. 10)  and  combine  (II. 9)  and  (II. 10) 


to  obtain 


H  =  l_i 

3o  3t  .2 
3y 


(IV. 35) 


where  <ji  is  the  plastic  strain  rate,  as  in  (II.  12a).  Substitution  of 
(IV. 34)  in  (IV. 35)  gives 


p  li  =  (S  A  +  S  9  +  SWP  ) 

3t  1  ,yy  2  ,yy  3  ,yy 

3S  3S  3S 

+  (  6  +  -  9  +  — -  WF  ) 

3y  ,y  3y  >Y  3y  >Y 


(IV. 36) 


where  the  coefficients  are  for  the  model  IV- 34,  given  by 


S_  =  — —  =  T/[niJ)  •  f  (n)  ] 

1  ~  *p 

3y 


(IV. 37a) 


S2  E  !?  =  -  aT  *  (n) ]/(202n2) 


(IV. 37b) 


e  _  3x  ,  1 

S_  =  -  =  mx  (  - - 1) , 

3  tt 


(IV. 37c) 


Equations  (II. 11)  and  (IV. 36)  constitute  the  governing  equations  for 
the  evolution  of  the  temperature  9  and  the  plastic  strain  rate  <f>  =  yP, 
x  in  these  equations  is  obtained  from  (IV. 34)  and  WP  is  obtained  from 
the  integration  of  the  second  part  of  (II. 5). 

For  a  homogeneous  deformation  the  coefficients  are  y-independent , 


S.  =  S. (t) 
l  l 


for  i  =  1,  2,  3 


(IV.  38) 


Moreover,  if  we  confine  attention  to  the  late  time  response;  then  we 
can  set 


s3<t  ->  »)  =  0 


(IV.  39) 


Based  on  these  assumptions  the  governing  equations  reduce  to 


p  =  s  <t>  +  s  e 

o  3t  1  ,yy  2  ,yy 


30 


*  =  r  0  +  r  t(0,<)>) 

3t  o  ,yy  1 


where 


t  (8  ,4>)  =  —  [f  (n)  ]C(6) 
/3 


Introduction  of  the  solution  vector 


HU 


allows  equations  (IV. 40. a, b)  to  be  written  in  the  form 


Z  =  A(2)  Z  +  A(0)  Z 

~.t  ~ ,yy  ~ 


where  the  coefficient  matrices  A^  and  A^°^  are  given  by 


(2) 


VP0  Yp0 


(0) 


'0 


r^t  0 


Taking  the  Fourier  transform  (cf.  (IV. 9))  of  equation 
obtains 


Z  =  c(t)  Z 

9  ^ 


where 


2  2 
-?S1  "5S2 


c(t)  = 


V  ro 


*  .  »  ,  •  »  ^  t  *  »  '  »  '  •  '  h'  *  •  m  *  \  '  »,  ■ 


(IV. 40. a) 

( IV. 40 .b) 

(IV.40.C) 

(IV. 41) 

(IV. 42) 

(IV. 43) 

(IV. 42)  one 

(IV. 44) 

(IV. 44. a) 
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For  the  stability  criterion  (IV. 22),  stability  requirements  for  the 
considered  model  can  be  given  by 
e2 

(|sj  +|S2|)  <  a  ;  a  >  0  (IV. 45. a) 

o 

(r  x  +  S2r  )  <  3  ;  6  >  0  ( I V . 4 5 . b ) 

1  o 

The  inequality  (IV.45.b)  is  always  satisfied;  however  the  inequality 
(IV. 45. a)  is  violated  as  soon  as  either 

yP  -*  2Dq  ,  i.e.  n  -►  1  (IV. 46. a) 

or 

6  ->--2-  ,  i.e.  n(0)  -*■  0  (IV.46.b) 

lb! 

The  second  limit  is  reached  for  all  homogeneous  deformations  without 
heat  loss,  if  these  deformations  are  carried  to  sufficiently  large 
strains.  Thus,  the  material  model  characterized  by  (IV. 34)  is 
especially  sensitive  to  temperature  increases.  For  the  numerical 
values  given  in  [13]  for  a  and  b  the  limit  (IV.46.b)  is  reached  at  a 
temperature  0  =  1.72;  for  adiabatic  deformation  from  room  temperature 
this  temperature  is  reached  at  a  strain  of  approximately  200%.  This 
limitation  on  homogeneous  deformations  for  the  model  (IV. 34)  will 
be  discussed  further  in  chapter  (V) . 
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V.  FINITE-DIFFERENCE  SOLUTION  OF  THE  NONLINEAR  PROBLEM 

In  order  to  obtain  a  better  understanding  of  the  onset  of  locali¬ 
zation  in  simple  shear  deformation,  we  seek  a  finite-difference 
numerical  solution  to  the  system  (11.9,12)  for  the  general  form  of  the 
constitutive  equation  (11.12).  Numerical  solutions  of  this  problem 
have  been  presented  to  interpret  torsional  Kolsky  bar  experiments 
[13 J ,  [15].  An  important  feature  of  these  solutions  is  the  intro¬ 
duction  of  a  geometric  imperfection  in  the  specimen.  The 
numerical  calculations  by  Costin  [12]  show  evidence  of  localized 
deformation  in  a  low-carbon  cold-rolled  steel  (CRS)  at  high  strain 
rates,  and  stable  deformation  for  a  hot-rolled  steel  (HRS)  of  the  same 
composition.  These  results  are  in  agreement  with  experimental  observa¬ 
tions.  Merzer  [13]  discusses  the  effect  of  thermal  diffusion  on  band 
spreading  in  Kolsky  bar  experiments.  He  shows  that  for  a  thermal, 
visco-plastic  model  of  an  aluminum  alloy  the  band  does  not  spread 
significantly  at  strain  rates  of  500  s  1  and  higher. 

The  objective  of  the  present  study  is  to  examine,  through  numerical 
solutions,  the  response  of  different  materials  to  different  types  of 
initial  perturbations  in  the  field  quantities,  for  various  combinations 
of  boundary  and  loading  conditions .  A  number  of  constitutive  equations 
is  considered  in  order  to  study  the  relative  effects  of  strain-hardening 
and  thermal  softening.  We  consider  the  following  system  of  equations: 


3v  _  1_  3t_ 

at  ”  p0  3y 

3yP  _  3v_ 

3t  3y 


(V.l) 


(V.2) 


.V 
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l.I.  i  a  ' 


."In  V' 


(V.3) 


f  ; 
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Elastic  effects  are  not  considered  in  this  model,  in  order  to  facilitate 
the  calculations;  however,  this  approximation  is  not  expected  to  affect 
the  results  significantly  since  elastic  strain  rate  is  much  smaller  than 
the  plastic  strain  rate.  Equation  (V.3)  has  a  quasi-linear  term,  while 
equation  (V.4)  is,  in  general,  fully  nonlinear. 

Before  presenting  the  finite  difference  formulation,  we  point  out 
that  we  seek  the  qualitative  behavior  of  the  solution,  which  makes  the 
stability  of  the  difference  method  our  primary  concern  and  the  method's 
accuracy  of  secondary  importance.  However,  accuracy  has  to  be  kept 
within  reasonable  limits  to  satisfy  the  consistency  requirement.  An 
essential  objective  in  this  formulation  is  to  be  able  to  resolve  the 
time  scale  on  which  major  changes  in  the  solution  occur.  This  makes  an 
explicit  formulation  preferable  to  an  implicit  one  since  small  time 
steps  are  required  for  good  time  resolution  and  the  implicit  schemes 
have  the  associated  difficulty  of  requiring  the  solution  of  a  system 
of  simultaneous  equations  at  each  time  step.  Thus,  the  unconditional 
stability  feature  of  implicit  schemes  is  less  attractive  for  the  type 
of  problems  considered.  Use  of  an  explicit  scheme  requires  the  choice 
of  a  "proper"  time  increment,  where  a  "proper"  time  increment  is  the  one 
that  maintains  scheme  stability  and  is  large  enough  to  give  good 
computational  e  f  f ic iency . 

The  difference  formulation  is  presented  for  the  system  (V.l-4) 
with  the  following  auxiliary  conditions. 
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A.  Initial  Conditions 

For  a  constant  overall  strain-rate,  initial  conditions  for  velocity, 
strain  and  strain-rate  are : 

v(y,o)  =  y,  yP(yfo)  =  yq i  ?P(y,o)  =  1. 

The  initial  condition  for  the  temperature  is  assumed  to  be 

6 (y ,0)  =  1  +  6p(y) 

where  6p(y)  is  a  small  perturbation.  The  perturbation  is  taken  to  be 

either  a  smooth  periodic  fluctuation,  or  a  step  perturbation  over  a 

narrow  region  near  the  center  of  the  slab  thickness. 

P  -P 

The  initial  conditions  on  Y  ,Y  and  0  are  used  to  calculate 
the  initial  stress  t(y,0)  from  equation  (V.4). 

Alternatively,  for  quasi-static  deformations  with  a  geometric 
imperfection,  we  take  the  thickness  in  the  direction  z  to  be  a 
function  T  (y) .  From  equilibrium  the  shear  stress  must  satisfy 

T (y,0)T (y)  =  Const.  (V.6) 

The  initial  temperature  is  taken  to  be  0(y,O)  =  1.  The  corresponding 

•  P  P 

initial  strain  rate  y  (y,Q)  and  initial  strain  y  (y,0)  for  which 

•  P  P 

(V.4)  is  satisfied  is  obtained  by  assuming  y  (y,0)  and  y  (y,0)  to 
be  proportional ,  i . e . 

P  *P 

y  (y/0)  =  yqy  (y,o) • 

Substituting  this  expression  and  (V.6)  into  (V.4),  and  requiring  the 
mean  strain  rate  to  be  unity,  we  obtain  the  constant  in  (V.6)  and  the 
initial  strain  rate  distribution. 
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B.  Boundary  Conditions 

We  can  prescribe  either  velocities  or  stresses  at  the  boundaries 
y  =  0  and  y  =  1.  For  most  of  this  analysis  we  assume  velocity 
boundary  conditions  of  the  form 


V(0,t)  =  0 

V(l,t)  =  1 


(V.  7) 


We  assume  that  the  deformation  occurs  without  heat  transfer  from  the 
slab.  This  condition  is  ensured  by  the  boundary  conditions 


u.u-o. 


(V.8) 


V.l.  Finite  difference  formulation 

Our  basic  formulation  is  concerned  with  the  dynamic  problem  in 
which  inertial  effects  are  retained.  We  confine  our  attention  in  this 
part  to  the  case  of  velocity  boundary  conditions. 

The  key  idea  in  this  analysis  is  to  rewrite  the  system  of 
equations  (V.l, 4)  in  a  modified  form  that  reflects  its  "parabolic" 
nature.  For  this,  differentiate  (V.4)  with  respect  to  y  to  obtain 


9t 

By 


(V.  9) 


where 


dip 

(V.lOa) 

3YP 

_ 

Be 

(V. 10b) 

■'W.v.n. y*’>y -V-v-y 

<  •  .i.A--  . .  -  k. 
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(V. 10c) 


From  the  compatability  relation  (V.2),  we  get 


•P  2 

3y  =  3_  ,3v  =  3_v 

3y  “  3y  l3y  3y2 


(V.ll) 


Substitution  of  (V.ll)  and  (V.9)  in  (V.l)  gives 


_  1  ,  £>  30  Hi, 

3t  '  V  1  3y2  2  3  *  K 


(V.12) 


This  form  of  the  momentum  equation  can  be  viewed  as  a  parabolic  equa¬ 
tion  for  the  velocity  v(y,t)  which  is  similar  to  the  well  known 
diffusion  equation 


3u  _  3  u 

3t  .2 
3x 


(V.13) 


but  with  variable  coefficient  S^(y,t)  and  lower  order  terms.  Sub¬ 
stitution  of  the  shear  stress  of  (V.4)  into  CV.3)  allows  the  governing 
equations  to  be  written  as  the  following  system  of  three  evolutionary 
equations . 


=  1_{S  lil  +  S  +  s  ^>L-} 

9t  p0  1  3y2  2  3*  3 


(V . 14) 


111  =  ii 

3t  3y 


(V.15) 


^-T  +  r,,|M0»YP»yP)  *  !?- 


(V.16) 


Introducing  the  solution  vector 


we  can  write  the  system  (V.14-V.16)  as 


(V. 17) 


y  yy 

0)  =  A*!*)  +  AJ  0) 

~,t  -,y  ~ .yy 


(V. 18) 


where 


(V .  19) 


(V. 20) 


We  introduce  the  following  terminology  for  finite  difference 
quantities 

h  =  Ay  ...  spatial  increment 
k  =  At  ...  time  increment 


x_,  =  x(jh,nk)  •••  Value  of  the  variable  x  at 

.  , th  .  _  th 

the  j  node  and  the  n  time 

level  (assuming  that  x(y,t) 

coincides  -  at  least  locally  - 

with  the  described  values  x° 


at  the  respective  nodes) . 


We  make  use  of  the  difference  operators 


D+(h)  ....  forward  difference  operator 
D  (h)  ....  backward  difference  operator 


DQ (h)  ....  central  difference  operator 


where 


D  (h)xn 
+  3 

D  (h)  x11 


Vh)xi 


(xj+rXj)/h 

(xn-xn  ) /h 
3  3-1 

-,-x"  >/2h 


x  -X, 
3+1  j- 


(V . 21) 


(V. 22) 


(V . 23) 


The  difference  operator  D+(h)D  (h)  =  D  (h)D+(h)  represents  the 
symmetric  second  order  difference  operator,  i.e. 


D+(h)D_(h)x^  =  (Xj+1-2x”+x^_1)/h2. 


(V. 24) 


We  will  use  the  operators  DQ (h)  and  D+(h)D  (h)  to  approximate 
first  and  second  order  spatial  derivatives,  respectively.  Also, 


D+(k)  will  be  used  to  approximate  the  time  derivatives. 


We  will  solve  the  evolution  equations  (V.14,16)  in  the  given  order 


using  the  updated  values  to  modify  the  calculation.  Therefore, 
difference  equations  can  be  given,  for  interior  nodal  points,  as 


follows : 


D  (k)vn 
+  D 


=  —  {S?  D  (h)D  (h)Vn  +  s"  Do(h)9n  +  .D  (h)  (yV)  (V.25) 


lj  + 


2  j  0 


3  j  0 


D+(k)  (yP)"  =  DQ(h)v"+1 


(V. 26) 


find 


D+(k)9^  =  rQD+(h)D_(h)e^  +  r1  lip  (6^  ,Yj+1 , Yj+1)  ]DQ (h)  v^+1. 


(V. 27) 
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The  shear  stress  t .  is  evaluated  from 

1 


n+1  .  ,  .n+1  ,  P  n+1  .  ...  n+1  . 

-  tHQj  ,  (Y  ) j  ,  lD0(h>Vj  ]>• 


(V.28) 


Spatial  derivatives  at  the  boundaries  are  approximated  by  a  second 
order  discretization  derived  by  the  method  of  undetermined  coefficients 
The  results  are 


3f 
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y=0 
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(V. 29) 
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“  <f  fN  -  2Vl  +  I  *  °(h2> 


(V .  30) 


92f 
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2fl  +  f2)/h  +  °{h) 


(V. 31) 


y=0 


32f 


“  +  0(h) 


(V. 32) 
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Since  the  boundary  conditions  are  prescribed  for  the  velocity  v(y,t) 
and  the  temperature  0(y,t),  we  do  not  need  to  use  (V.3l,32).  There¬ 
fore  the  order  of  accuracy  of  the  scheme  (V.25,30)  becomes 

€  =  0{h2,k)  i.e.  0(2,1).  (V.33) 

Consistency  of  the  scheme  with  the  original  system  (V.18)  requires 
that  €  is  at  least  0(1,1),  therefore  (V.33)  implies  consistency  as 
well. 


V.2.  Stability  of  the  difference  scheme 

The  notion  of  scheme  stability  we  use  is  a  modification  of  the 
Von-Neumann  stability  criterion  [l 8 ]  discussed  by  Strang  [20]  which, 
for  the  system  (V.18) ,  gives  the  stability  requirement: 


k 


h 


2 


Max  eigenvalue  of 


1 

2  ' 


From  (V.20)  this  requirement  is 

(y,t) 


k  (Vy< 
(JL)maJ 

h  (  P0 


In  general,  S^^t)  >  k/c  so  that  the  limit  on  the  time  step  is 


k  < 


poh 


-  2S1(y,t) 


(V. 34) 


The  criterion  (V.34)  used  to  bound  the  time  increment  k  is  an 
approximate  one  based  on  freezing  the  coefficient  S^(y,t)  between 
time  levels.  We  evaluate  the  right  hand  side  of  the  inequality  at  each 


time  level  and  set 
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kn+^  =  (0.97) *min 

3  j 


poh2/2 


(Sl)j 


(V. 35) 


n  +  1 

where  k_.  is  the  time  increment  between  levels  (n)  and  (n+1)  . 


V. 3.  Formulation  for  stress  boundary  conditions 

The  analysis  of  the  case  where  stresses  are  prescribed  at  the 
boundaries  is  simplified  by  rewriting  the  governing  equations  in  a 
form  that  allows  the  solution  to  be  obtained  following  the  same  procedures 
used  in  the  case  of  velocity  boundary  conditions.  Differentiating  the 
constitutive  equation  (V.4)  with  respect  to  time,  one  obtains 


il.  s  ifll  +  s  +  s  ill 

at  bi  2  b2  at  b3  at 


(V.  36) 


where,  S1,S2  and  S3  are  defined  by  (V.10,a,b,c) .  From  (V.l)  and 


(V.2)  we  have 


3 JL_  1  32t 
3t2  p0  3y2 


(V .  37) 


Then,  substitution  of  (V.3),  (V.2),  and  (V.37)  into  (V.36)  gives  the 
following  evolutionary  equation  for  the  flow  stress 


3t  S1  32T  .  326  T  9V\  r-  3v 

3t  P0  3y2  2  °  3y2  1  3y  3  ** 


(V .  38) 


Introducing  the  solution  vector 
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(V . 39) 
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we  can  write  the  governing  equations  in  the  form 


z  «  BYz  +  B^z 

~,t  ~,y 


(V. 40) 


where 


S  /p 

r  0 


VS2  +  S3 


(V .41) 


(V.42) 


The  numerical  solution  of  equation  (V.40)  is  obtained  by  an  ex¬ 
plicit  finite  difference  scheme  analogous  to  that  described  in  the 
previous  section.  The  spatial  derivatives  are  replaced  by  central 
differences  while  time  derivatives  are  replaced  by  forward  differences. 
The  system  of  evolutionary  equations  is  solved  in  the  order  of  the 
components  of  the  vector. 

The  scheme  stability  is  analogously  determined  by  the  requirement 


Max  eigenvalue  of  B 
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(V. 43) 
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Computational  evidence  show  that  scheme  stability  is  ensured  by  using 

»oh2 

k  =  0.97  (— — )  (V . 44) 

2S1 

which  is  the  same  step  size  used  for  the  case  of  velocity  boundary 
conditions . 

V.4.  Formulation  for  the  case  of  geometric  inhomogeneity 

We  seek  a  finite-difference  solution  to  the  general  nonlinear 
system  (11.9,12)  for  the  case  of  a  thin-walled  tube  with  a  groove  in 
it  as  shown  in  Fig.  (V.l).  The  tube  is  subjected  to  torsion  which,  for 
sufficiently  small  values  of  the  ratio  of  wall  thickness  to  tube  radius, 
is  a  good  approximation  of  simple  shear. 

Numerical  simulation  of  this  problem  has  been  presented  by 
Litonski  [151,  Costin,  et  al,  [12]  and  Merzer  [13]  for  the  case  of 
quasi-static  deformation.  They  assumed  homogeneous  deformation  con¬ 
ditions  in  each  uniform  thickness  section  of  the  tube.  Merzer  [13] 
included  heat  conduction. 

A  complete  numerical  solution  for  this  problem  is  presented  — 
accounting  for  inertial  effects,  and  without  imposing  the  assumption 
of  deformation  homogeneity  in  uniform  thickness  regions.  The  system  of 
governing  equations  is  taken  to  be  (V.l, 4).  Due  to  the  geometric 
discontinuity,  modifications  are  required  in  the  initial  conditions. 

In  addition,  boundary  conditions  are  required  at  the  surface  of  dis¬ 
continuity.  To  obtain  these  conditions,  consider  a  strip  of  width  26 
at  the  edge  of  the  groove  as  shown  in  Fig.  (V.2). 

Let  and  '+'  denote  evaluations  of  the  various  quantities  at 
the  left  and  right  hand  sides  of  the  edge  shown  in  Fig.  (V.2). 


H 
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Figure  (V.l)  -  Longitudinal  Cross-section  in 
the  thin-walled  tube  with  the  groove  (6,  S) 


Figure  (V.2)  -  Imaginary  strip  at  the  border  of  the  groove 
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Continuity  of  particle  velocity,  temperature,  and  shear  force 
(or  torque)  gives 


v  =  v 


0=0 


T  =  At  ;  A  <  1 


(V.45a) 
(V.45b) 
(V. 45c) 


where  A  =  A  /A  is  the  ratio  of  the  cross-sectional  areas. 

Consideration  of  energy  flux  into  the  strip  of  width  2£,  and 
taking  the  limit  as  6  -*■  0  gives  the  additional  relations 


(!V  =  wiV 

3y  M3y'  * 


(V. 45d) 


Equations  (V.45a-d)  are  the  conditions  that  are  to  hold  at  the  interface 

y  ■  y«* 


V.4.a.  Finite-difference  solution 

The  symmetry  of  the  problem  allows  us  to  consider  only  the  solution 
in  the  half  y  €  [0,1/2].  The  initial  conditions  are  chosen  such  that 
the  system  of  governing  equations  (V.l-4)  and  the  continuity  require¬ 
ments  (V.45)  are  met.  Plastic  strain  and  temperature  are  taken  to  be 

.P 

homogeneous.  Plastic  strain  rate  y  is  assumed  to  jump  at  the  inter¬ 
face  in  order  to  give  the  required  jump  in  the  flow  stress.  The 
overall  strain  rate  is  maintained  constant. 

The  difference  method  is  the  same  as  described  in  section  V.l, 
except  for  the  treatment  of  the  geometric  discontinuity  at  y  =  yQ. 
Discretization  at  nodes  j+  and  j  (see  Fig.  (V.3))  is  performed 
by  replacing  the  spatial  derivative  3/3y  by  forward  and  backward 


\-Wv\ 

■V.y.VA.V. 


differences  analogous  to  (V.29)  and  (V.30).  The  velocity  at  level 


n  +  1  is  obtained  by  proceeding  from  the  boundary  y  =  0  at  which 
v(0,t)  =  0  to  the  nodal  point  j  -  1  using  equation  (V.25)  and 
similarly  proceeding  from  the  boundary  y  =  y  at  which  v(j,t)  =  j 
to  the  nodal  point  j  +  1.  In  order  to  evaluate  (v?+1)  we  use  the 
jump  relations  (V. -,5b)  and  (V.45c)  and  the  constitutive  equation  (V.4) 


x-  *(0  »YP.f  In-)  -  *ce  ,yp  ,  |~) 


j  o 


+  '  dy 


j-  8y 


(V. 46) 


where  all  quantities  are  to  be  evaluated  at  the  level  n  +  1.  The 
velocity  gradients  in  (V.46)  are  approximated  by  the  one-sided  ex¬ 
pressions  (V.29),  (V.30).  With  the  continuity  of  the  velocity  field 
at  y  =  y  equation  (V.46)  becomes  an  equation  for  the  velocity  at 

^  I  ^ 

the  interface  v^  .  This  equation  is  solved  by  either  taking  0 ^ 

P  P 

and  y  ,Y  to  be  the  values  at  the  level  n  or  by  an  iterative 
j+  3" 

P  P 

procedure  in  which  the  values  of  6 . , y  . ,Y  are  updated  as  they  are 

D  j" 

obtained  from  (V.26)  and  (V.27),  with  centered  differences  replaced  by 


one-sided  differences. 


The  temperature  is  evaluated  by  means  of  equation  (V.27)  for  all 
nodal  points  except  the  nodes  j+,j  at  the  interface.  The  temperature 
at  these  nodes  is  obtained  from  (V.45b)  and  (V.45d) ,  using  one-sided 
difference  operators.  Finally,  the  flow  stress  is  obtained  from 
equation  (V.28). 


V.5.  Numerical  results 

The  numerical  solution  to  the  fully  nonlinear  system  (V.1-V.4)  is 
presented  for  a  number  of  constitutive  equations  modelling  the  response 


of  visco-plastic  materials.  Initial  imperfections  of  the  homogeneous 
solution  are  taken  as  a  periodic  temperature  perturbation  or  as  a 
permanent  geometric  defect  (cf.  V.4).  Numerical  results  are  presented 
for  both  velocity  as  well  as  traction  boundary  conditions.  The  dy¬ 
namic  problem  is  considered  for  all  models  discussed.  An  essential 
objective  is  to  compare  numerical  results  with  linear  analysis 
predictions  (Chapter  IV)  in  order  to  determine  the  usefulness 
of  these  predictions. 


A.  Velocity  boundary  conditions 

We  present  herein  numerical  results  for  the  case  where  velocities 


are  prescribed  at  the  boundaries  as  given  by  conditions  (V. 7) . 


A.l.a  Power  law 

Consider  the  constitutive  equation 

x  =  (0)V(YP)m(-yp)n  (V.47) 

where 

v  =  thermal  softening  (hardening)  exponent 
m  =  strain  hardening  exponent 
n  =  strain-rate  exponent 

A  series  of  numerical  solutions  has  been  obtained  for  n=l  and  for  different 
values  of  v  and  m  in  order  to  study  the  relative  effects  of  thermal 
softening  and  strain-hardening.  Most  visco-plastic  materials  exhibit 
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strain-hardening  and  thermal  softening.  Under  such  circumstances. 


it  is  believed  that  at  high  strain  rates  deformation  instability  occurs 


Hi 


as  soon  as  the  net  hardening  rate  vanishes,  where  'A'  refers  to 

dt*A 


an  adiabatic  deformation. 

From  the  numerical  solutions  it  appears  that  the  vanishing  of  the 
adiabatic  hardening  rate  is  not  sufficient  to  ensure  the  onset  of 
localization.  Indeed,  the  numerical  solution  for  n=l  shows  remarkably  stabl^ 
deformation  even  when  the  slope  of  the  isothermal  stress-strain  curve 
becomes  negative . 


Numerical  results  have  been  obtained  for  the  following  values  of 
the  model  parameters  corresponding  to  an  aluminum  alloy. 


Thermal  properties 


k 

po 
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-1  -1 
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The  overall  strain  rates  considered  are  in  the  range 


,3  -1  6  -1, 

[10  s  -  10  s  ] . 


The  initial  homogeneous  strain  is  taken  as  yQ  =  0.01,  while  the 


perturbation  in  the  initial  temperature  distribution  is  chosen  as 


yy) 


€  sin [2ir (—  -  y) ] 


(V.48) 


where  €  =  0.005  or  G  =  0.05.  The  dimensional  value  of  the  homo¬ 


geneous  (unperturbed)  initial  temperature  distribution  is  0Q  =  300°K. 


The  slab  thickness  is  H  =  0.01m.  The  homogeneous  initial  value  of 


8 


the  flow  stress  is  taken  as  t_  -10  Pa. 
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The  dimensionless  spatial  step  size  is  chosen  as  Ay  =  1/40,  while 
the  temporal  step  size  is  monitored  throughout  the  calculation  to  meet 
the  stability  requirement  (V.35). 

Numerical  solutions  for  the  power  law  (V.47)  with  m  >  0,  v  <  0,  and 
n  =  1  show  a  slight  departure  from  the  homogeneous  solution  at  small 

p 

strains  y  <0.05.  Then,  the  deviation  from  the  homogeneous  solution 
reaches  an  upper  bound  at  which  the  inhomogeneous  solution  retracts  back 
towards  the  homogeneous  solution,  but  at  much  slower  rate.  Figure  (V.4) . 

For  the  material  considered,  the  rate  of  temperature  rise  is  slow,  with 
only  a  1  percent  temperature  rise  for  a  5  percent  plastic  strain  incre¬ 
ment.  The  stabilizing  effect  of  thermal  diffusion  seems  to  be  almost 
negligible  in  the  deformation  considered,  Figure  (V.6) .  An  increase  in 
the  amplitude  of  the  initial  inhomogeneity  to  6  *  0.05  causes  a  pro¬ 
portionate  increase  in  the  excursion  from  the  homogeneous  solution,  but 
it  does  not  change  the  qualitative  features  of  the  solution,  i.e.  an 
envelope  is  attained  and  the  deformation  remains  stable. 

Finite-difference  solutions  have  also  been  obtained  for  cases  with 
-1  <  m  <  0  and  v  <  0  for  which  elementary  linear  stability  analysis  (e.q.[ll]) 
allows  the  possibility  of  exponential  growth.  For  cases  of  thermal 
softening  and  strain  softening  with  n  =  1,  the  flow  stress  decreases 
monotonically  as  the  solution  evolves  vrhile  the  numerical  solution  continues 
to  become  more  inhomogeneous  as  time  increases ;  however ,  the  time  scale 
on  which  significant  departure  from  the  homogeneous  solution  occurs  is 
extremely  long — apparently  too  long  for  the  localization  to  correspond 
to  the  shear  band  development  observed  in  experiments.  A  typical  situ¬ 
ation  is  illustrated  in  Figure  (V.5)  which  corresponds  to  v  =  -0.75, 
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m  =  -0.95  and  n  =  1  at  a  dynamic  global  strain  rate  of  (10  )s  .  Clearly, 
the  profiles  at  later  times  are  getting  closer  to  each  other.  This 
appears  to  be  due  to  a  decrease  in  the  plastic  working  as  a  result  of 
the  rate  of  decreasing  stress  values.  It  appears  that  the  solutions  for 
(m  <  0)  tend  asymptotically  to  the  limiting  situation  of  a  localized 
deformation.  However,  if  we  confine  ourselves  to  strain  amplitudes  that 
are  commonly  obtained  in  the  plastic  deformation  of  real  materials,  we 
reach  the  conclusion  that  the  deformation  described  by  the  model  (V-47) 
is  remarkably  stable  for  n  =  1. 

The  accuracy  of  the  finite  difference  code  has  been  examined  by 
calculating  a  solution  for  a  special  case  (v  =  m  =  0,  n=l)  for  which 
an  exact  solution  is  available.  Close  agreement  with  the  exact  solution 
was  obtained. 


A.l.b  Arrhenius  law 

Consider  the  constitutive  equation 

i 

T  =  tb(Y  )  -  <— )*n(— ) 

y 


(V. 49) 


discussed  in  Chapter  III. 

Numerical  results  have  been  obtained  for  the  same  Aluminum  alloy 
as  in  A.  Other  constants  required  in  (V.49)  are 


a)  =  10 


12  -1 
s 


^  =  1.38  x  10"23J/°K 
W  =  0.5  x  10_26m3 


( V.  50) 


while  the  strain  hardening  function  has  the  form 


tb(YP)  =  tb<0)  •  El  +  b(l  -  e"CY  )) 


(V.  51) 


where 


108  Pa 
0.5 


c  =  4.0. 

We  also  introduce  the  following  additional  dimensionless  quantities 
(cf.  [14] ) 

A  A  P 
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(V. 52) 
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Numerical  results  obtained  for  this  deformation  are  shown  in 
Fig.  (V. 8) .  The  qualitative  nature  of  the  deformation  behavior  is 
similar  to  model  (A. 1. a)  for  m  <  0  and  v  <  0,  in  the  sense  that  the 
solution  has  a  tendency  to  become  more  inhomogeneous  as  time  increases 
however,  the  deviation  from  the  homogeneous  solution  again  occurs 
very  slowly.  The  deformation  governed  by  the  Arrhenius  law  remains 
contained  even  when  the  point  of  maximum  stress  is  reached,  which 
appears  to  be  strong  evidence  for  stability.  In  this  sense  the  numer¬ 
ical  results  appear  to  be  consistent  with  the  predictions  of  the 
linear  stability  analysis  by  Clifton  et  al.  [14]. 

The  numerical  result  shown  in  Fig.  (V. 8)  is  obtained  using  the 
hardening  function  given  by 


& 
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p  ”  *T)  P  P  * 

T0<Y  >  =  <  "—  *  >  *  Y  +  fQ  ,  for  Y  S?  Y 

Y 


(V. 53) 


.  P  * 
,  for  Y  >  Y 


where  the  dimensional  values  of  the  parameters  are 


fQ  =  10 8  pa  ,  f  *  1.5  x  lo8  pa 


(V.  54) 


Y  =0.15 


A.l.c  Bodner-Merzer  law 


Consider  the  constitutive  equation 


T. (WP) 

t  =  — -  [f (n) ]c(6> 

/I 


(V. 55) 


defined  in  (IV. 5. 3).  Numerical  values  for  the  material  constants  are 


taken  to  be  the  same  values  given  by  Merzer  [13].  Non- uniformity  of 


the  strain-rate  along  the  length  of  the  specimen  continues  to  increase 


as  the  deformation  proceeds.  For  e  =  0.05  the  deformation  becomes 


strongly  localized.  Computed  profiles  for  e  =  0.005  are  shown  in 


Figs.  (V.9.a-d). 


The  instability  obtained  for  the  model  (V. 55)  can  be  attributed  to 


two  main  reasons : 


1)  The  strong  thermal  softening  that  overcomes  strain  hardening 


throughout  the  deformation  process,  see  Fig.  (V.9.c); 


2)  The  extreme  sensitivity  to  temperature  which  has  been  dis¬ 


cussed  in  (IV. 5. 3) . 


The  instability  reported  by  Merzer  [13]  for  this  model  -  in  the 


case  of  a  geometric  imperfection  -  has  been  ascribed  to  the  "persistency" 


of  the  geometric  defect;  however  in  view  of  the  unstable  response 
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of  this  model  to  a  "transient"  temperature  perturbation  it  is 
evident  that  the  details  of  the  constitutive  model  as  well  as-the 
nature  and  magnitude  of  the  initial  imperfection  play  important  roles 
in  the  localization  process. 


7 
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A. 2  Geometric  imperfection 

We  consider  numerical  results  for  the  problem  described  in  (IV. 4) 
in  which  a  permanent  groove  is  introduced  in  a  thin-walled  tube,  Fig. 

(V. 1) .  The  groove  is  characterized  by 

A  =  1.04  and  6  =  0.1. 

A. 2. a  Power  law 

The  computed  profiles  for  the  model  (V. 47)  are  shown  in  Figs. 

(V.lO.a-d)  and  (V. ll.a-d).  Clearly,  for  n  =  1  the  power  law  -  again  -  gives 
strongly  stable  response  in  spite  of  the  initial  jump  in  the  strain  rate, 
Fig.  (V.lO.d) .  The  difference  between  temperature  levels  inside  and 
outside  the  groove  remains  less  than  one  percent  up  to  total  strains 
of  40  percent.  Figure  (V.lO.b). 

For  the  case  of  m  >  0;  the  stress  rises  monotonically  while 
it  decreases  monotonically  for  m  <  0,  Figures  (V.lO.c),  (V.ll.c)  re¬ 
spectively.  The  maximum  rise  of  the  strain  rate  from  its  initial  value 
(inside  the  groove)  is  less  than  one  percent  which  is  strong  evidence 
of  the  stability  obtained  for  the  power  law  given  by  (V.47) ,  with  n  =  1. 
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The  effect  of  strain-rate  hardening  on  the  response  for  the  power 
law  given  by  (V.47)  is  examined  by  choosing  the  parameters  v,  m  and 
n  to  correspond  to  two  steels,  namely  a  cold-rolled  steel  (CRS-1018) 
and  a  hot-rolled  steel  (HRS-1020) .  Experimental  investigations  [19] 
showed  that  deformation  localizes  into  shear  bands  in  the  CRS  whereas 


the  HRS  showed  a  stable  response. 

Numerical  values  for  the  parameters  used  in  this  computation  are 


CRS-1018 


HRS-1020 


V 

=  -0.38 

V 

=  -0.51 

m 

=  +0.015 

m 

=  +0.12 

n 

=  +0.019 

n 

=  +0.0133 

T 

=  436  MPa 

T 

=  261  MPa 
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Thermal  Properties: 
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The  length  of  the  specimen  is  taken  as  2.5  mm  and  a  dynamic 
global  strain  rate  of  10^  s  ^  is  considered.  The  groove  characteristics 
are 

X  =  1.02  and  6  =  0.05 

The  results  of  this  calculation  are  shown  in  Figures  (V.12.a-f)  and 
(V.13.a-f)  for  the  HRS  and  the  CRS  respectively.  The  solutions  indicate 
remarkably  distinct  behavior  for  the  two  steels  considered.  The 
strong  hardening  rate  exhibited  by  the  HRS  seems  to  stabilize  the 
deformation,  whereas  the  deformation  of  the  CRS  localizes  inside  the 
groove  as  soon  as  thermal  softening  overcomes  strain  and  strain-rate 
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hardening.  Furthermore,  for  the  CRS  the  shape  of  the  profiles  in¬ 
side  the  groove  indicates  a  strong  localization  process.  These  re¬ 
sults  are  -  qualitatively  -  in  agreement  with  experimental  observa¬ 
tions  in  [19] .  The  instability  observed  for  the  CRS  with  n  <<  1  and 
the  stability  shown  in  Figures  (V.ll)  for  n  =  1  are  consistent  with 


>. 
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where  "L.O.D."  stands  for  "lower  order  derivatives".  Differentiating 
(V.58)  w.r.t.  time,  we  get  the  following  wave  equation  for  the  plastic 
strain  rate 


2  2 
3  4  r  1  a  A 

rf  ■  ‘r  (s3 +  w1  ijf +  L-°-D- 

d  t  O 


(V. 59) 


Observing  that  <  0  for  the  thermal  softening  features  of  the 


Arrhenius  law,  we  can  see  that  at  early  times  when  can  be  neg¬ 


lected  with  respect  to  S^,  a  wave- like  response  is  obtained  when  the 


strain  hardening  (S^  >  0 )  dominates  the  thermal  softening  (S^  <  0)  . 


This  analysis  is  supported  by  the  numerical  results  shown  in  Fig. 

(V. 16. a-d)  where  no  hardening  is  considered  and  the  wave-like  early 
response  is  totally  eliminated. 

Consideration  of  the  late  time  response  for  the  case  of  the 
Arrhenius  law,  suggests  that  the  existence  of  a  permanent  groove  in 
the  specimen  contributes  to  an  instability.  This  can  be  seen  by 
noting  the  large  temperature  increase  inside  the  groove  relative  to 
that  outside  it,  as  well  as  the  increase  of  strain  rates  inside  the 
groove  and  their  decrease  outside  the  groove.  In  spite  of  these 
instability  features  we  regard  this  response  as  essentially  stable 
for  practical  purposes.  That  is,  although  the  deformation  begins 


with  strain  rates  inside  the  groove  which  are  4.5  times  larger  than 

P 


those  outside,  the  final  strain  rates  at  y  =1.0  are  only  6  times 
larger  than  those  outside  the  groove.  Therefore  the  observed  in¬ 
stability  is  very  weak  even  with  a  substantial  initial  jump  and  with 
continuation  of  the  deformation  to  large  strains.  We  prefer  to  reserve 
the  term  instability  for  cases  analogous  to  reported  experimental 


results  on  shear  bands  in  which  the  strains  within  the  bands  are  much 
larger  than  in  neighboring  regions. 

The  results  shown  in  Figure  (V.15.a-d)  are  obtained  for  a 
shallower  groove,  X  =  1.004,  the  difference  between  temperature  levels 
in  and  out  of  the  groove  is  less  pronounced.  Moreover  the  difference 
between  strain  rates  is  less  than  that  obtained  for  X  =  1.04,  and 

p 

the  strain  rate  inside  the  groove  at  y  =1.0  remains  well  below  the 

P 

value  at  y  =0.0. 

The  results  obtained  for  the  non-hardening  idealization  are  pre¬ 
sented  in  Figure  (V.16.a-d) .  Wave-like  early  response  is  completely 
eliminated  upon  removal  of  the  strain  hardening.  At  late  times  the 
differences  in  strain  rate  are  comparable  to  those  obtained  when 
strain  hardening  was  included.  This  result  is  consistent  with 
regarding  the  response  as  stable  for  the  case  of  the  Arrhenius  law, 
since  even  for  such  strongly  destabilizing  conditions  as  a  groove  and 
no  strain-hardening  the  deformation  does  not  become  much  more  non- 
uniform  at  late  times  than  it  is  initially.  Computational  experiments 
performed  for  the  Arrhenius  law  revealed  negligible  sensitivity  of 
the  results  to  heat  conduction. 

A.2.c  Bodner-Merzer  model 

Numerical  results  for  the  model  (V.55)  have  been  reported  by 
Merzer  [13]  for  a  specimen  with  a  groove.  Our  results  for  the  same 
problem  are  shown  in  Figs.  (V.17.a-d).  In  order  to  improve  the 
accuracy  of  the  solution  near  the  groove  a  finer  mesh  is  used  than  in 
[13]  and  jumps  across  the  groove  boundaries  are  accounted  for  ex- 


The  behavior  for  m  >  0  is  remarkably  stable.  The  deformation 
tends  to  become  more  homogeneous  with  increasing  time,  see  Fig.  (IV. 18. a) 
This  stable  response  results  even  when  the  point  of  maximum  load  is 

P 

passed,  i.e.  for  y  >0.025.  On  the  other  hand,  the  response  for 
m  <  0  is  remarkably  unstable.  Shear  bands  form  at  the  two  boundaries 
where  strain  rates  are  24  times  higher  than  their  initial  values  (see 
Fig.  (V.19.d)).  Moreover  this  intense  shearing  occurs  at  relatively 

p 

small  strains  (y  =  0.05) ,  indicating  a  strong  bifurcation  process. 

As  a  result  of  thermal  softening  and  strain  softening,  the  flow 
stress  decreases  monotonically  except  at  y  =  0  and  y  =  1.  Sustaining 
high  stress  levels  at  the  boundaries  leads  to  the  softening  of  the 
material  there  and  to  the  resulting  shear  strain  localization. 

A  similar  qualitative  response  is  observed  for  the  Arrhenius  law 
given  by  (V.49)  with  no  hardening.  The  deformation  tends  to  become 
more  inhomogeneous  as  time  grows;  however,  the  speed  of  deformation  is 
slower  than  that  obtained  for  the  power  law,  see  Figures  (V.20.a-d). 
Strain  rates  are  only  (1.5)  times  their  initial  values  at  50%  of  plastic 
straining.  This  slow  localization  can  be  explained  by  the  weak  thermal 
softening  exhibited  by  the  Arrhenius  law  as  well  as  the  absence  of  the 
isothermal  strain-softening,  see  Figure  (V.20.c). 

The  response  for  the  Bodner-Merzer  constitutive  model  is  -  again  - 
qualitatively  similar  to  earlier  results;  however,  the  rate  of  heat 
generation  is  much  higher  as  compared  to  the  power  law  or  the  Arrhenius 
law,  see  Figure  (V.21.b) .  The  limit  (IV.46.b)  imposed  on  temperature 
by  the  constitutive  law  (V. 35)  is  reached  at  the  specimen  boundaries  for 
total  plastic  strain  of  approximately  (35)  percent. 
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The  remarkable  feature  of  the  solution  of  the  simple  shear  problem 
with  traction-controlled  boundaries  is  that  the  same  resnonse  is  obtained 
whether  or  not  the  initial  solution  is  perturbed.  In  other  words, 
the  solution  of  an  unperturbed  dynamic  problem  with  stress  boundary 
conditions  becomes  immediately  inhomogeneous  as  deformation  proceeds. 

The  evolution  of  inhomogeneous  deformation  seems  to  be  a  consequence 
of  the  imposed  stresses  at  the  boundaries,  since  holding  stresses  to 
be  constant  there  while  allowing  stress  to  decrease  in  the  interior 
creates  sharp  stress  gradients  in  the  vicinity  of  specimen  boundaries; 
as  a  result,  strain  rates  become  large  -  forming  regions  of  localized 
deformation.  This  feature  of  the  problem  with  stress  boundary  condi¬ 
tions  is  discussed  further  in  Appendix  (III) . 

V. 6  Concluding  remarks 

For  velocity  boundary  conditions,  simple  shearing  deformations 
are  remarkably  stable  for  the  three  constitutive  equations  considered 
except  for  the  Bodner-Merzer  model  and  the  power  law  with  n  <<  1. 

The  localization  observed  for  the  Bodner-Merzer  model  appears  to  be 
due  to  a  non-physical  degeneracy  of  the  model  at  finite  temperatures 
that  may  be  reached  during  adiabatic  experiments.  The  localization  ob¬ 
served  for  the  CBS  with  n  <<  1  is  viewed  as  indicative  of  the  effects 
of  thermal  softening,  weak  strain  hardening  and  weak  strain  rate 
sensitivity  in  obtaining  significant  strain  localization  at  moderate 
strains.  Although  a  persistent  geometric  defect,  such  as  a  groove, 
appears  to  be  slightly  more  destabilizing  than  an  initial  perturba¬ 
tion  in,  for  example,  temperature,  the  essential  conclusions 
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regarding  the  stability  of  simple  shearing  deformations  do 
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not  appear  to  be  changed  substantially  by  the  type  of  inhomogeneity. 


However,  the  strength  of  the  inhomogeneity  plays  a  significant  role. 

For  stress  boundary  conditions  it  appears  that  simple  shearing 
deformations  are  much  less  stable.  For  a  power  law  model  with  iso¬ 
thermal  strain  softening,  localized  deformations  occur  at  relatively 
small  strains. 
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Figure  (V.9.a)  -  Initial  temperature  perturbation,  name  parameters 
as  in  Merzer's  Calculation,  c«0.005,  profiles  arc  obtained  for 
intervals  of  (yp»0.05). 


Figure  (V.O.b)  "  Temperature  profiles  for  the  problem 
defined  in  rig.  (V.9.a). 
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Figure  (V.9.C)  -  stress  profiles  for  the  nrohlem 
defined  in  rig.  (V.9.a). 
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Figure  (V.9.d)  -  strain  rate  protiles  for  the  problem 
defined  in  Fig.  (V.h.a). 
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Figure  (V.lO.a)  -  Geometric  imperfection,  6*0. 1,  A-1.04,  v— 0.75, 
nm+0.25,  £-0.0,  plotted  values  are  identical  in  the  range  yp*0.0-0.4 
for  the  shown  scale. 


Figure  (V.lO.b)  -  Temperature  profiles  for  the  probl 
defined  in  Fig.  (V.lO.a). 
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Figure  (V.lO.c)  -  stress  profiles  for  the  problem 
defined  in  Fig.  (V.lO.a). 


Figure  (V.lO.d)  -  strain  rate  profiles  for  the  problem 
defined  in  Fig.  (V.lO.a). 
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Fiqure  (V.12.a)  -  Geometric  imperfection,  rt=0.05,  >*1.02, 
v—0.51,  m=0.12,  n=0.013,  plots  are  obtained  for  intervals 
of  yP=0.05. 


Fiqure  (V. 12 . b)  -  Temperature  profiles  for  the  proble 
defined  in  Fig.  (V.12.a). 
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Fiqure  (V. 12.d)  -  Strain-rate  profiles  for  the 
problem  defined  in  Fiq.  (V.l2.a). 
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Figure  (V.13.a)  -  Geometric  imperfection,  6=0.05,  A«1.02, 
v*-0.38,  m*0.015,  n=0.019,  plots  are  obtained  for  intervals 
of  yPsr0.05 


Figure  (V.13.b)  -  Temperature  profiles  for  the 
problem  defined  in  Fig.  (V.13.a). 
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Figure  <V.13.c)  -  Stress  profiles  for  the  problem 
defined  in  Fig.  (V,13.a),  stress  decreases  after  y  =0.1. 


Figure  (V.ll.d)  -  strain-rate  profiles  for  the  problem 
defined  in  Fig.  (v.l).a),  beyond  yf>=0.05?  strain-rate 
increases  rapidly  inside  the  groove  whereas  decreasing 
outside  of  it. 
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Figure  (V.H.e)  -  Plastic  strain  profiles  for  the 
problem  defined  in  Fig.  (V.13.a). 


CRS-1018 


0.0?  0,12  0, 16 
average  strain 


0,20  0.211 


Figure  (V.ll.f)  -  Computed  Nominal  stress-strain  relation j  stress  attains 
its  maximum  value  at  approximately  yP”0.1 
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Figure  (V.lS.a)  -  same  a»  rio.  (V.14.a)  except  for  X-1.004, 
profiles  obtained  for  yp>0.2  tend  to  move  towards  the 
initial  profile  (yp«0.0). 


^  '  ^Hiperacure  profiles  for  the  probl 
defined  in  Fig .(V.lS.a),  profiles  are  obtained  for 
intervals  of  (?P-0.05). 
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Figure  (v.ls.c)  -  stress  profiles  Tor  the 
problem  defined  in  Fig.  (V.lS.a). 


Figure  (V. 15. d)  -  strain-rate  profiles  for  the  problem 
defined  in  Fig.  (v.lS.a)i  a  similar  response  to  that 
in  Fig.  (V.14.d)  is  observed. 
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Figure  (V.lS.a)  -  Initial  temperature  perturbation, 
c-0.005,  strees  boundary  conditions,  v«-0.75,  m*0.1, 
n»l;  profiles  are  obtained  for  intervals  of  ^*0.005. 


Fiaure  (V.lS.b)  -  Temperature  profiles  for  the 
problem  defined  in  Fin.  (V.lS.a). 


POWER  LOW,  (STRESS  B.C.) 


TP’0.025 


rP-o. 005 


nj.oo  o.zo  o.vo  o.eo  o. Bo  1.00  i.zo 

SPRTIflL  DIMENSION 


Flour*  (V.lB.c)  -  stress  profile*  for  the  problem 
defined  In  Fig.  (V.lS.a)  -  stress  reaches  Its 
maximum  at  (fp«0.  025) . 
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Figure  (V.ls.d)  -  strain  rate  profiles  for 
the  problem  defined  in  l-'in.  (V.lS.a). 
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Figure  <v.21.a)  -  Unperturbed  initial  solution,  stress 
boundary  conditions,  plots  are  obtained  for  intervals 
of  yP-0.05. 


Figure  <V.21.b)  -  Temperature  profiles  for  the  problem 
defined  in  Fig.  (V.21.a),  notice  that  the  temperature 
values  at  the  boundaries  reach  the  limit  (IV.4$.b) 
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Figure  <v.21.c)  -  Stress  profiles  for  the  problem 
defined  in  Fig.  (v.21.a> . 


Figure  (V.21.d)  -  strain-rate  profiles  for  the  problem 
defined  in  Fig.  (V.21.a). 


APPENDIX  (I) 


Closed  Form  Solution  of  a  Linear,  1  Order 


Time-Dependent  System  of  Equations 


Consider  the  system: 


F  =  W(t)F 


with  the  initial  condition 


?<v  -  ? 


Claim: 


If  W(t)  commutes  with  W(x)dx, 


(  ft 

F(t)  =  exp  j  W(x)dx  F(tQ) 


is  the  unique  solution  of  (1)  and  (2) 


Proof : 


Let  K(t)  be  a  fundamental  matrix  of  the  system  (1) ,  i.e. 

K(t)  =  W ( t)  K(t)  (4 

and  det.  K(t)  /  0.  Then  the  solution  of  (1)  and  (2)  can  be  written 


F(t)  =  K (t)  K_1 (tQ)  F° . 


Introduce  the  notation 

♦<tft0)  =  K(t)  K_1 ( tQ)  I 

for  the  fundamental  matrix  <j>(t,tQ)  which  is  the  solution  operator 


of  the  system  (1)  and  (2) . 


Assume 


Where 


4> (t , tQ)  *  exp  |  M(t,tQ)  | 
M(t, t„)  =  [  W ( x) dx 


It  is  clear  that  4> ( t , tQ)  satisfies  (1),  i.e. 

$(t,t  )  =  W(t)  4>(t,t0) 

We  seek  a  proof  that  (7)  satisfies  (9) .  Moreover,  we  are 
that  this  solution  is  unique,  (Haar's  uniqueness  theorem)* 

exp{M(t)  }  =  I+  M  +  ^-MM+...  . 

Then  the  R.H.S.  of  equation  (9)  becomes 

R.H.S.  =  W  +  WM  +  —■  WMM  +  •••  . 

While  the  L.H.S.  of  (9)  becomes 

L.H.S.  =  rjy  t4>  (t,tQ)  ] 

=  [I  +  M  +  -r-—  MM  +  •••] 
at  i ! 

2 

=  M  +  yy  MM  +  • • • 

*  W  +  ~  MW  +  ~  M^W  +  •  •  • 

=  W  +  MW  +  yy  MMW  +  . 

Comparing  expressions  (10),  (11);  we  conclude  that,  if  w(t) 
with  M(t) ,  i.e. 

MW  =  WM 

then  (10)  and  (11)  are  identical,  and  consequently  4> (t , tQ ) 


(9) 

guaranteed 


(10) 


(11) 

commutes 


is  Jie 


unique  solution  to  (9). 


>  V, 
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APPENDIX  II 


Numerical  Solution  of  Merzer's  Problem 


The  outlines  of  the  numerical  algorithm  are  discussed  in  (V.4). 
The  governing  equations  used  are 


.  _  «5  3ZV  .  e  30.  +  c  3yP 

P0  3t  “  S1  .  2  S2  3y  S3  3y 
dy 

3yP  3v^ 

3t  “  3y 

36  ,  K  t  326  .  6  _  3yP 

w  * 


T1(WP) 


ltC-l/2)/(a/0+b)] 


2  An(“i 
Y 


where 


where 


S  s*I - 1 _ 

1  3yp  nYP.f(n) 


S„  =  =  -ax-in[f(Ti)]/(2e2n2) 

A  0  V 


S,  =  ~~jp  m  at/u-a) 
3Y 


A  -  mx  (- - 1) 

1 


,  P,  -  .  -mW 

T1(W  5  =  Ti  "  <VTo,e 

x .  lf(nn=w> 
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f(n)  *  2  in  n  ;  n 

c(8>  ■  sw  ■  nl9> 


substitution  of  (6)  and  (7)  in  (5)  gives  a  nonlinear  equation  for  v 
which  is  solved  by  Newton's  method.  The  initial  velocity  distribution 
is  obtained  in  a  similar  manner  where  the  jump  in  the  flow  stress  is 
satisfied  through  an  equivalent  jump  in  the  strain  rate  assuming  con¬ 
tinuous  strain  and  temperature  distributions. 
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Finite-difference  formulation 

Computational  experiments  for  the  considered  system  of  equations, 
with  Merzer's  values  for  the  parameters,  indicate  that  the  formal 
explicity  scheme  is  no  longer  efficient.  Therefore,  the  scheme  is 
improved  by  using  the  two-level,  explicit,  stable  Du-Fort  and  Frankel 
algorithms,  with  careful  choice  of  k,  h  to  maintain  consistency.  The 
additional  initial  value  is  evaluated  by  the  formal  explicit,  one- 


level  method;  the  formal  method  with  time  step  k1  is  used  until  we 
reach  t  =  nk.  =  k_;  then  the  calculation  is  continued  by  the  Du-Fort 


and  Frankel  method  with  a  uniform  time  step  k2  >>  k^ .  The  efficiency 
of  the  algorithm  is  strongly  improved  by  this  technique.  The  accuracy 
is  also  improved  since  the  Du-Fort  method  is  second  order  in  time 


and  space. 


The  Du-Fort  and  Frankel  approximation  to  the  standard  heat  equation 


3u  _  3  u 

3t  =  V  3x2 


is  given  by 


n+1  n-1  n  ,  n+1  n-1.  n 

ui  -  U1  Vi  ■  lui  4  ui  1  *  Vi 


Our  system  has  variable  coefficients  and  lower  order  terms;  however, 
computation  evidence  shows  the  utility  of  this  algorithm.  The  algorithm 
described  in  (V.4)  is  applied  for  this  system  of  equations  with  the 
exception  that  we  use  Merzer's  constitutive  equation  and  the  differ¬ 
encing  procedure  illustrated  in  (9) . 


V*  "  •  *  O  s.  •».  *.  O  *.  s  v  A  *.  *•  *•  *•  *•  -   -  /•  •  .  •  .  •-  *  .  .  %  *  -  *  .  •  - 
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APPENDIX  III 


On  the  lack  of  Homogeneous  Solutions  for  the  Case  of 
Stress  Boundary  Conditions 

Consider  the  momentum  and  compatability  equations  (II. 9)  and  (II. 10) 


\n 


3y  9v  1  3t 

3t  3y  G  3t  1  ' 

Assume  there  exists  a  homogeneous  solution  for  the  velocity  and  stress 
given  by 


v(y,  t)  =  V(t) y 
x(y,  t)  =  T ( t)  . 

Substituting  solution  (3) ,  (4)  in  (1)  ,  (2)  we  obtain 


pQ  V(t)y  =  0  -*■  V(t)  =  0 


V(t)  =  V  (Constant) 
o 


Moreover,  Equation  (2)  yields 


yP  -  V  -  T(t)  I 

o  g  'velocity  b.c. 


Now,  let  us  distinguish  two  cases 


A.  Velocity  boundary  Conditions: 


In  this  case,  V  (0 ,  t)  =0  ,  v(l,  t)  =  Vq  ,  and  Equation  (6)  is 


•P  i  A  , 

=  V  -  ;=■  T(t)  I 

o  G  'velocity  b.c. 


(7. A) 


where  the  last  term  is  the  stress  rate  obtained  for  the  homogeneous 


solution  corresponding  to  the  velocity  boundary  conditions. 


B.  Stress  Boundary  Conditions: 

In  this  case,  T(t)  is  prescribed  at  the  boundaries  and  -  in  turn  - 
across  the  specimen  since  the  solution  is  homogeneous,  hence;  Eguation  (6) 
reads 


Yp  =  Vo  -  i  ;<t) 


prescribed 


(7.B) 


Comparison  of  Equations  (7. A)  and  (7.B)  allows  us  to  conclude  that  the 
only  way  by  which  (7.B)  is  satisfied  becomes 

T  ( t)  |  ..  =  T  ( t)  | 

'prescribed  velocity  b.c. 

i 

In  other  words,  the  unique  possible  homogeneous  solution  is  the  one 
that  corresponds  to  velocity  boundary  conditions,  otherwise  the  solution 

becomes  inhomogeneous  as  deformation  proceeds. 


Norm  in  the  Linear  Stability  Analysis 


Consider  the  system  of  linearized  partial  differential  equations 


defined  in  (IV. 7)  by 


L(o))  =  0 


This  system  can  be  rewritten  as 


A  oj,.  =  B_U3,  +  B,o),  +  B_o) 

~  t  2~  yy  1~  y  o~ 


where 


B2  =  -AYy,  Bx  =  -Ay  and  BQ  =  -A 


Introduce  the  change  of  variables 


where 


oj  =  Dz 


Then  the  system  of  equations  (2)  becomes 
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2’t  -  B2-'yy  +  Cl?'y  *  V 


where 


=  DB^D  and  CQ  =  DBQD 


The  matrices  B2  and  are  symmetric;  B2  is  positive,  semi -definite. 

Now  consider  the  rate  of  change  of  the  energy  norm  of  the  solution 


fH 

!?»?>  =  ?•;« 

J  n 


From  (4)  and 


we  obtain 


^<z.z>  -  <(B2z,yy  ♦  C1?.y  +  C0Z),Z>  * 

+  (b2?'w  *  ci?'v  *  co5l'> 


Integration  of  the  second  derivatives  by  parts  gives 


<B25'yy'?>  +  <5'B25'yy^=  ~2^B2?'y'5'y^  (8) 

plus  boundary  terms  which  vanish  for  boundary  conditions  of  either  pre¬ 
scribed  heat  flow  or  prescribed  temperature.  Since  B2  is  positive, 
semi-definite,  the  right  side  of  (8)  is  non-positive.  Integration 
of  the  first  derivatives  in  (7)  by  parts  gives 


^C^z,  ,z^  +  (z,C^z,  y  =  -^C^z,z,y)  +  (z,CjZ,^ 


plus  boundary  terms  which  vanish  for  prescribed  velocities  or  prescribed 


shear  tractions.  From  the  symmetry  of  C^,  the  right  side  of  (9)  is  zero 


Use  of  the  results  (8)  and  (9)  in  (7)  gives 


^<z,z>  <  2<Cqz,z> 


(10) 


or 


(ID 


where  |j  z|  =  ^z,z^  and  1 1 (t )  1 1  2  is  the  L2~norm  of  the  matrix  CQ(t) 


Integration  of  (11)  gives 


|z(t) | |2  <  { exp I 


[0nc° 


(n) I l2dn 


2(0) 


(12) 


or 


ls(t)||2  <  |exp|^  max  ^  |  | CQ (n)  |  l2tj|  j  |z  (0)  |  |. 


(13) 


Using  the  equivalence  of  norms  for  finite-dimensional  spaces,  we  can,  for 
convenience,  replace  the  L2~norm  of  Cg  with  the  maximum  norm  to  obtain 

|  | z  (t)  |  1 2  <  jexp|a  max|  |cQ(n)  |z(0)  |  |2  (14) 

The  well-posedness  of  the  system  (2)  is  ensured  provided  that 

max  | |c  (n) | I  <  o  (15) 

ne[0,t]  ° 

where  o  is  a  positive  constant. 

The  estimate  (15)  gives  the  inequalities  (IV. 24)  obtained  from  the 
analysis  in  Chapter  IV  for  finite  wave  numbers  K • 


